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I.  INTRODUCTION: 


This  ’’jport  presents  the  clutter  information  obtained  during 
the  clutter  measurement  task  performed  in  support  of  the  SEEK 
IGLOO^Radar  Development  Program.  The  SEEK  IGLOO  Program  is  di- 
rected toward  the  implementation  of  minimally-attended  radar 
systems  at  the  Alaskan  Air  Command  surveillance  radar  sites.  The 
successful  implementation  of  unattended  radars  would  significantly 
reduce  the  high  cost  of  maintaining  personnel  at  each  site.  How- 
ever, such  systecis  require  automatic  detection  schemes  with  low 
false  alarm  rates  and  moderate  probabilities  of  detection  in  all 
clutter  environments.  Therefore,  accurate  knowledge  of  the  clutter 
environment  is  necessary  to  obtain  a properly  designed  radar. 

The  primary  purpose  of  the  clutter  measurement  task  was  to 
determine  the  worst  case  amplitude  statistics  and  spectra  of  land 
clutter.  Radar  backscatter  was  measured  from  range-azimuth  loca- 
tions characterized  by  three  terrain  types:  valleys  (muskeg),  hills 
and  mountains.  Both  high-sited  and  low-sited  radars  were  used  in 
the  measurement.  Other  measurements  performed  include  sea  clutter 
amplitude  and  spectra,  weather  clutter  spectra,  and  angel  spectra. 

This  report  begins  with  a background  discussion  on  the  measure- 
ment team  composition  and  the  general  procedure  at  each  site.  The 
measurement  technique  and  equipment  is  discussed  in  Section  III 
followed  by  the  presentation  and  discussion  of  the  measured  ampli- 
tude clutter  in  Section  IV  and  spectral  data  in  Section  V. 
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II. 


BACKGROUND: 


The  measureiuent  team  was  comprised  of  individuals  from  the 
Electronic  Systems  Division  (ESD)  , 4754th  Radar  Evaluation  Squad- 
ron (RADES) , Rome  Air  Development  Center  (RADC) , and  Alaskan  Air 
Command  CAAC) . The  functions  and  members  of  each  organization 
are  given  in  Figure  1. 

The  clutter  data  was  collected  at  five  (5)  selected  AAC 
radar  sites  during  the  period  2 August  through  17  September  1976. 
These  sites  were  selected  for  their  severe  or  unusual  clutter 
environments.  Table  I gives  the  location  of  the  sites  and  the 
reason  for  selection. 

The  schedule  at  each  site  is  given  in  Figure  2.  During 
the  first  three  days,  the  FADES  and  AAC  teams  calibrated  the  radar. 
This  procedure  was  performed  to  assure  properly  operating  equip- 
ment prior  to  the  measurements  and  to  provide  a check  of  the 
critical  antenna  electrical  tilt  and  other  system  parameters  for 
use  in  the  measurement  analysis.  These  calibrations  typically 
.ncluded  the  following: 

Antenna  System 
Mechanical  Checks 

Solar  Measurements  (except  for  Tatalina,  Cape  Newenhamj 

Transmitter  System 

PRF 

Pulse  Shaping 
Power  Output 
Frequency 


FIGURE  1:  Team  Composition 
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FIGURE  2:  Schedule 


Receiver  (Normal  and  MTI) 

Bandwldths 
Center  Frequencies 
MTI  System 
Video  Levels 

Plan  Position  Indicator  (PPl) 

Focus  and  Stability 
Azimuth  Orientation 
Range  Mark  Accuracy 

In  addition,  the  RF  losses  between  the  receiver/transmitter  to 
the  antenna  feed  and  the  mixer-preamplifier  dynamic  range  were 
measured. 

Various  measurements  were  performed  during  the  last  four 
days.  These  included  PPI  clutter  photos,  amplitude  cumulative 
distribution  measurements,  and  spectral  measurements.  Due  to 
operational  considerations,  nearly  all  measurements  were  performed 
between  1800  - 0600  local  time  and  the  operational  downtime 
allowed  for  spectral  measurements  was  typically  9 to  12  hours. 
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m.  MEASUREMENT  EQUIPMENT  AND  METHODOLOGY: 


A.  Amplitude  Measurements; 

The  clutter  parameters  measured  were  the  effective  radar 
cross-section  o^  and  the  effective  coefficient  of  reflection  as 
defined  in  Table  II.  The  term  "effective"  is  used  because  the 
measured  parameter  contains  an  unknown  propagation  factor  F which 
was  not  factored  out.  The  equipment  used  in  obtaining  the  ampli- 
tude cumulative  distribution  function  of  these  parameters  was  an 
implementation  of  the  radar  range  equation.  Briefly,  the  equip- 
ment had  the  following  features: 

1.  R’**  (or  R"*)  range  normalizing  thresholds  were  available. 

2.  Op  Cot  o')  reference  could  be  related  to  the  threshold 
using  the  radar  parameters. 

3.  Op  (or  Oq)  amplitude  profiles  could  be  obtained  by 
inserting  attenuation. 

4.  The  range -azimuth  measurement  interval  was  programmable. 

5.  Threshold  crossings  were  sampled  at  an  interval  equal 
to  or  greater  than  once  per  pulse  width. 

6.  A sample  pulse  width  of  100  NS  was  used  to  increase 
independence  between  adjacent  samples. 

7.  Threshold  crossings  were  counted  using  an  electronic 
counter. 

8.  A threshold  output  for  PPI  display  was  available  for 
pictures . 
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TABLE  II;  Definition  of  Oq  and  Oq 


rt  p**  m 

PrR-(4fl)*  is 
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PtGtGRl* 

^ - 

<Jc  F-  . r PrR*I*sC4Tr)*'| 

r /r 
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Effective  Radar  cross-section 

Radar  cross-section 

Propagation  factor 

Received  power  at  input  to  receiver 

Transmitter  peak  power 

Range 

Gain  of  transmitting  antenna 
Gain  of  receiving  antenna 
Wavelength 

System  RF  loss  - LjLrLpLra 
Effective  coefficient  of  reflection 
Coefficient  of  reflection 

Resolution  area  * R ( — ) SinQs  Sec<j)g  ( ) 

2 /T 

Speed  of  propagation 
Pulse  width 

Antenna  beamwidth  (One  Way) 

Transmitter  RF  loss 
Receiver  RF  loss 
Propagation  loss 
Radome  loss 
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9.  One  measurement  per  antenna  scan  of  one  area  could  be 
taken. 

A generalized  block  diagram  of  the  clutter  amplitude  analyzer 
(CAA)  is  shovn  in  Figure  3.  A block  diagram  of  the  threshold  gen- 
erator and  comparator  is  given  in  Figure  4.  The  RF  log  amplifier 
had  a dynamic  range  of  over  80  dB  with  a log  accuracy  within  ^ 1 
dB.  The  DC  log  amplifier  had  a dynamic  range  over  80  dB  with  a 
maximum  log  conformity  error  of  1%.  The  linear  ramp  had  a peak 
deviation  from  linear  of  less  than  5%  over  the  range  used.  By 
varying  the  slope  a,  gain  y,  DC  bias  B,  and  time  delay  Z,  the 
threshold  was  adjusted  to  give  the  desired  threshold-versus- time 
curve  such  as  those  shown  in  Figure  5.  The  reference  parameters 
®c  ref  °o  ref  related  to  the  thresholds  by  the  equations 

derived  in  Appendix  A. 

The  methodology  used  was  straightforward.  Areas  of  strong 
clutter  were  located  on  the  PPI  and  correlated  with  terrain  type 
using  topographic  maps.  Regions  with  range  and  azimuth  extents 
within  the  terrain  type  were  outlined  and  converted  into  the  pro- 
grammable parameters  of  the  measurement  equipment  CFigure  6). 
Received  signals  were  amplified  and  compared  to  the  threshold. 
Since  the  threshold  represents  a constant  clutter  parameter  value 
for  the  given  radar  parameters,  those  signals  exceeding  the  thresh 
old  represent  regions  with  a higher  clutter  parameter  value.  The 
threshold  comparator  output  was  then  sampled  at  a rate  of  once  per 
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FPS-83A 


FIGURE  3:  Block  Diagram  of  CAA  Interface  with  FPS-93A 


THRESHOLD 

COMPARATOR  ► THRESHOLD  OUTPUT 


FIGURE  4:  Threshold  Generator  and  Comparator 




! V"'"/ 


X H$ 


V— 


Nomogmaoiis  Tirrain 


y > 

G f Azimuthal  distance  between  samples  ■ — ®P 

A - Start  Range  ■ 

B = Stop  Range  = Ngxx 

R ■ B - A “ (Ng  - Nyy)  Tjj  = Range  Extent  of  Samples 
Tx  ® Distance  in  range  between  samples  i (.1  usee)  ^ t 

e * Total  azimuth  extent  of  samples  = NO 


CNb  - Na)N 


Total  number  of  samples  taken 


R,  PRF,  e,  T,  Op  are  known  constants  for  a given  measurement 
and  N0,  NAt  Nb»  Nt»  programmable  parameters.  The  sampling 

interval  is  started  manually  or  automatically  at  azimuth  x. 
FIGURE  6:  Example  of  CAA  Operation 
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independent  range -azimuth  cell  for  the  programmed  number  of 
samples  corresponding  to  the  range-azimuth  extent  of  the  area 
being  measured.  The  sampled  outputs  that  exceed  the  threshold 
were  counted.  Then  the  probability  that  a radar  return  from  the 
measured  area  will  be  greater  than  the  threshold  was  estimated  as 

p 'i  /v,  Total  No.  of  Threshold  Crossings 

r (.a  > _ 

Total  No.  of  Samples 

By  taking  such  measurements  with  increments  of  attenuation  in- 
serted, thereby  changing  a cumulative  distribution  function 

was  easily  plotted  for  each  area  of  interest. 

The  typical  pulse  widths  used  in  the  measurements  were  0.5, 

1,  3 and  6 microseconds.  All  amplitude  measurements  were  taken 
using  an  IF  bandwidth  of  approximately  Z MHz  to  provide  information 
usable  for  certain  pulse  compression  systems  (i.e.,  13  bit  Barker 
coded  6 usee  pulse)  as  well  as  a bounding  value  for  non-pulse  com- 
pression 0.5  - 6 usee  designs.  To  give  more  statistical  signifi- 
cance to  the  measurements,  ten  measurements  of  each  measured 
sector  were  performed  at  each  5 dB  attenuation  level  and  averaged 
to  obtain  each  plotted  value.  All  measurements  were  made  with 
threshold- to- rms  noise  ratios  of  10  dB  or  more  within  the  measure- 
ment range  to  reduce  the  contamination  caused  by  noise.  Finally, 
the  average  antenna  gain,  which  varies  as  a function  of  terrain 
and  distance,  was  calculated  for  each  measurement  to  reduce  the 
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error,  A detailed  error  analysis  of  the  measurements  is  presented 
in  Appendix  B. 

B.  Spectral  Measurements: 

The  purpose  of  the  spectral  measurements  was  to  obtain 
information  on  the  intrinsic  radar  clutter  spectra  of  land,  sea, 
weather  and  if  possible,  angels.  The  spectrum  analysis  equipment 
consisted  of  a model  440  Mini-Ubiquitous  Portable  Spectrum  Analyzer 
(Nicolet  Scientific  Corp.)  and  radar  interface  equipment  as  shown 
in  Figure  7,  This  analyzer  performs  an  FFT  on  1024  time  samples 
and  displays  400  frequency  points  on  a display.  The  resolution 
and  maximtua  display  frequency  was  selectable  from  .025  Hz  at  10  Hz 
to  50  Hz  at  20  KHz,  Other  useful  features  of  the  analyzer  include 
4 averaging  modes  and  a choice  of  linear  or  log  displays.  The 
information  was  recorded  on  oscillographs  of  the  spectrum  analyzer 
output.  The  spectrum  analyzer  oscillograph  annotation  is  given 
in  Figure  8,  The  interface  was  designed  primarily  for  analyzing 
frequencies  less  than  50  Hz  ('v>12  knots)  but  could  be  extended  up 
to  one-half  of  the  PRF  by  accounting  for  the  low  pass  filter  and 
sample-and-hold  attenuation  given  in  Figure  9. 

The  methodology  used  was  straightforward.  Areas  were  chosen 
for  examination  thru  observation  of  the  MTl  PPI  display.  These 
areas  were  correlated  with  a topographic  map  to  determine  the 
terrain  type.  The  antenna  was  stopped  and  set  to  the  chosen 
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HETRQDYNE 


FIGURE  7:  Block  Diagram  of  Spectrum  Analyzer  Interface 


The  above  template  is  provided  to  decode  the  photographs 


that  follow: 


1. 

Amplitude  Controls  (dB) 

7. 

Amplitude  Scale: 

2. 

Max  Frequency  Displayed 

VLN  - Linear 

VLG  - Log 

3. 

Display  Mode; 

8. 

Amplitude  Display  Max  CdB) 

S ■ Averaging  Mode 

Cursor  Frequency 

I - Instantaneous  Mode 

T - Time  Display 

9. 

10. 

Hz 

4. 

Number  of  Spectra  Averaged 

11. 

Amplitude  of  Display  at 

5. 

Data  Weighting:  W - Hanning 

Cursor  Frequency 

U - Uniform 

12. 

dB 

6. 

Frequency  Scale: 

HLN  - Linear 

HLG  - Log 

FIGURE  8:  Spectrum  Analyzer  Annotation  Identification 
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FIGURE  9:  Cascaded  Sampler -And- Ho Id  and  LPF  Transfer  Function 


azimuth,  "searchlighting"  the  clutter  area  to  be  examined.  The 
received  signals  were  heterodynea  to  baseband,  amplified  and 
filtered.  A chosen  range  cell  within  the  linear  range  of  the 
receiver  was  boxcarred,  providing  an  audio  signal  for  the  spectrum 
analyzer. 

First,  the  degree  of  spectral  contamination  caused  by  trans- 
mitter noise  and  system  instability  was  bounded  by  examining  the 
spectra  of  several  "stable"  clutter  returns  from  fixed  point  tar- 
gets at  different  ranges.  Then  chosen  clutter  areas  at  the  same 
or  lesser  ranges  were  examined  and  recorded. 


IV. 


SPATIAL  AMPLITUDE  DISTRIBUTIONS  OF  TERRAIN  AND  SEA  CLUTTER: 


The  probability  distributions  for  land  clutter  are  shown  in 
Figures  10  through  32.  The  presentation  of  this  data  is  broken 
down  into  five  (S)  categories: 

a.  Flat  lands,  such  as  swamps,  muskeg,  plains,  Figures  10-11. 

b.  Hills  and  mountains,  pcrtially  wooded  - Figures  13  - 23. 

c.  Volcanic  mountains  - Figures  24  - 26. 

d.  The  composite  distribution  of  all  terrain  features  over 

a 360®  azimuth  from  15  to  50  and  15  to  100  nautical  miles 
range  - Figures  27  and  28. 

e.  Comparative  distributions  of  different  locations  of  flat 
lands,  hills,  and  mountains  - Figures  29  - 32. 

One  (1)  sea  clutter  distribution  is  presented  in  Figure  33. 

These  figures  show  the  percentage  of  the  samples  in  which  the 
reflectivity  exceeds  the  amplitude  in  dB  shown  on  the  abscissa. 

Also  shown  on  each  graph  is  the  site,  date,  time,  pulse  width 
and  a description  of  the  terrain  feature  which  the  data  represents. 
Except  when  noted,  each  figure  shows  a family  of  curves  with  pulse 
width  as  a parameter.  Since  each  site  is  currently  part  of  the 
Alaskan  early  warning  system,  area  codes  are  given  in  lieu  of  the 
range-azimuth  locations  of  each  measurement. 

All  the  figures  indicate  that  the  distribution  functions  range 
between  log  normal^ and  contaminated  normal^ Examples  of 
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these  distributions  are  shown  in  Figure  34  with  a Rayleigh  dis- 
tribution for  comparison  and  the  medians  of  each  normalized  to 
zero  dB. 

From  an  examination  of  the  data,  the  distributions  tend 
toward  log-normal  as  the  sample  areas  became  uniform  or  homoge- 
neous. This  is  brought  out  by  the  straight  line  nature  of  the 
distribution  for  muskeg  backscatter.  The  effect  is  shown  partic- 
ularly in  the  tails  of  Figures  10  through  13.  Also  note  that  for 
this  type  of  terrain  at  low  grazing  angles,  such  as  at  Fort  Yukon 
(Figure  10)  and  Cape  Newenham  (Figure  12) , the  log-normal  spread^ 
between  the  median  and  the  84  percentile  is  large,  from  11  to  16 
dB.  For  the  higher  grazing  angles,  such  as  at  Tatalina,  Figure 
11,  overlooking  a sloping  wooded  area  broadside,  the  spread  is 
lower,  about  7 dB.  This  is  perhaps  due  to  less  shadowing^. 

Note  the  spread  also  tends  to  be  larger  for  shorter  pulse  widths. 
The  spread  for  the  predicted  log-normal  distribution  was  scaled 
from  the  curves  by  extrapolating  the  bottom  (lower  percentile) 
end  toward  the  top  (SO  percentile)  end. 

For  hilly  terrain  with  average  elevation  increasing  with 
range,  the  curves  begin  to  deviate  from  a straight  line  and 

appear  as  contaminated-normal.  This  begins  to  occur  in  the 

/ 

curves  for  Cape  Newenham  (Figure  12),  Indian  Mountain  (Figure  13), 
and  at  Tatalina  (Figure  21).  All  three  of  these  areas  have  hills 
rising  in  the  direction  away  from  the  radar  and  have  relatively 
uniform  granularity. 


As  larger  areas  of  shadowing  and  attendant  facing  slopes  are 
observed,  a pronounced  knee  appears  in  the  curves.  This  is  brought 
out  by  the  sector  observed  at  Tatalina  (Figure  22)  where  approxi- 
mately 20%  of  the  area  is  shadowed  and  about  501  of  the  area  is 
sloped  such  that  the  radar  grazing  angle  is  large. 

Extreme  cases  of  contaminated-normal  are  seen  in  sectors 
observed  at  Indian  Mountain  (Figure  15)  where  deep  valleys  and 
rising  mountain  faces  extend  over  large  portions  of  the  sector. 

The  two  (2)  curves  in  this  figure  are  very  irregular  suggesting 
large  differences  in  the  contaminated-normal  parameters^. 

As  one  observes  high  mountains  over  very  large  sectors  such 
as  at  Indian  Mountain  (Figures  16  and  17)  and  at  Tatalina  (Figure 
23),  the  irregular  effects  in  the  curve  appear  to  decrease.  This 
can  be  due  to  the  large  areas  under  observation  and  the  averaging 
of  the  reflectivity. 

The  composite  curves.  Figures  27  and  28,  show  a very  large 
difference  in  slope  between  the  50  percentile  and  the  low  percen- 
tiles but  the  transition  is  smooth.  Here  the  averaging  effects 
described  in  the  preceding  paragraph  are  occuring  over  the  entire 
curve.  For  example,  the  averaging  of  the  individual  curves  in 
Figures  29  to  31  would  tend  toward  this  result. 

The  conclusions  which  can  be  drawn  from  the  data  and  discus- 
sion are  that: 
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1.  Homogeneous  clutter  whose  scatterers  are  larger  than  a 
carrier  wavelength  and  much  less  than  a resolution  cell 
has  a log-normal  (straight  line)  distributed  reflectivity. 
The  smaller  the  grazing  angle  the  larger  the  spread 
between  the  50  and  84  percentile. 

2.  The  more  heterogeneous  the  data  becomes,  where  large 
areas  within  a sector  become  shadowed  and  large  exposed 
areas  become  eminent,  the  more  contaminated-normal  like 
the  distribution  becomes  with  very  irregular  knees  and 
bends.  Note  we  are  citing  large  areas  of  shadowing  and 
exposure,  not  numerous  small  areas. 

3.  As  the  areas  of  shadowing  and  exposure  become  numerous 
by  observing  larger  and  larger  sectors,  the  distribution 
takes  on  a smooth  transition  from  a low  slope  nt  the  50 
percentile  to  a high  slope  at  the  low  percentiles. 

4.  For  the  worst  case  conditions  at  the  tails  of  the  dis- 
tribution curve  we  may  assume  a log-normal  function 
according  to  Table  III.  This  will  give  an  upper  bound 
for  low  false  alarm  thresholds.  An  approximate  per  cent 
occurrence  of  the  three  terrain  types,  obtained  from 
contour  maps  of  several  sites,  are  also  given. 
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TABLE  III:  Worst  Case  Land  Clutter  Model 
(Extrapolated  from  the  Tails  of  Distributions) 


PER 

EXTRAPOLATED 

REFLECTIVITY  DISTRIBUTION 
AT  TAILS  OF  CURVES 

TERRAIN 

TYPE 

VEGETATION 

HEIGHT 

(FT) 

CENT 

OCCUR- 

RENCE 

MEDIAN 

84TH 

PERCENTILE 

Mountains 

Bare 

4000 

13 

0 

1 

-14 

Wooded 

Grass-Moss 

40v/0 

53 

-19 

-12 

Hills 

Small  Trees 

1000 

-- 

— 

— 

Valleys/ 

Heavily  Wooded 

1000 

34 

-20 

-13 

Muskeg 

Spruce,  Birch 

— 

— 

— 

Everything 

(Composite: 

15-100  NM) 

-23.5 

-17.5 

Z7 


r ustc  TERRAIN  TYPE  AREA  NO. 


TATALINA 

TATALINA 


REFLECTIVITY  { OB ) 

FIGURE  11:  Swamp  (Tatalina) 


CURVE  NO. 


REFLECTIVITY  ( DB  ) 

FIGURE  13:  Hills  (Indian  Mountain) 


PLACE  • -DATE 
CAPENEWENNAM  4 SEPT 


TERRAIN  TYPE  AREA  NO 


REFLECTIVITY  ( DB  ) 

FIGURE  16:  Hills  (Indian 


TERRAIN  TYPE  AREA  NO 
MOUNTAINS  6 


REFLECTIVITY  ( 08  ) 

FIGURE  17:  Mountains  (Indian  Mtn) 


REFLECTIVITY  ( OB  ) 

FIGURE  20:  High  Hills,  Mountains  C^^ort  Yukon) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE  AREA  NO 

15  TATALINA  25  AUG  6.0  LOW  MTNS/HIGH  HILLS  9 

18  TATALINA  25  AUG  3.0  LOW  MTNS/HIGH  HILLS  9 

19  TATALINA  25  AUG  1.0  LOW  MTNS/HIGH  HILLS  9 


REFLECTIVITY  ( DB  ) 

FIGURE  21;  High  Hills,  Mountains  (Tatalina) 


T usee  TERRAIN  TYPE  AREA 


FIGURE  23 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE  AREA  NO. 

1 COLD  BAY  14  SEPT  5.7  VOLCANIC  MOUNTAINS  1 

3 COLO  BAY  14  SEPT  5.7  VOLCANIC  MOUNTAINS  t 

5 COLO  BAY  14  SEPT  3.0  VOLCANIC  MOUNTAINS  1 

8 COLD  BAY  14  SEPT  1.0  VOLCANIC  MOUNTAINS  1 


CURVE  NO.  PLACE  DATE  r usee  TERRAIN  TYPE  AREA  NO. 

9 COLD  BAY  15  SEPT  0.5  VOLCANIC  MOUNTAINS  1 


FIGURE  25:  Volcanic  Mountains  (Cold  Bay) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE  AREA  NO 

2 COLD  BAY  14  SEPT  5.7  VOLCANIC  MOUNTAINS  2 

4 COLD  BAY  14  SEPT  5.7  VOLCANIC  MOUNTAINS  2 

6 COLO  BAY  14  SEPT  3.0  VOLCANIC  MOUNTAINS  2 

7 COLO  BAY  14  SEPT  1.0  VOLCANIC  MOUNTAINS  2 


FIGURE  26:  Volcanic  Mountains  (Cold  Bay) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE  AREA  NO 

3 TATALINA  23  AUG  6 EVERYTHING  3 


REFUCTIVITY  ( OB  ) 

FIGURE  27:  Composite  Distributions  (Tatalina) 


CURVE 


REFLECTIVITY  ( DB  ) 

FIGURE  28:  Composite  Distributions  [Tatalina) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYI 

12  TATALINA  25  AUQ  6 HILLS 

13  TATALINA  25  AUG  5 SWAMPS.  RIVER 

14  TATALINA  25  AUG  6 McKINLEY  RA 

15  TATALINA  25  AUG  6 LOW  MTS/HIGH 


vssiosav  SNI0333X3  JO  AiniOVOOiM 
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REFLECTIVITY  ( DB  ) 

FIGURE  29:  Comparative  Distributions  (Tatalina) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE  AREA  NO 

16  TATALMA  25  AUG  3 RIVER  BEDS  & SWAMP  7 

17  TATAUNA  25  AUG  3 HILLS  6 

18  TATALiNA  25  AUG  3 MOUNTAINS  9 


REaECTIVITY  ( OB  ) 

FIGURE  30:  Comparative  Distributions  (Tatalina) 


CURVE  NO.  PLACE  DATE  T usee  TERRAIN  TYPE 

to  TATAUNA  25  AUG  1 MOUNTAHUS 

20  TATALINA  25  AUG  t HHiS 


REFUCTiVmr  ( DB ) 

FIGURE  31:  Comparative  Distributions  [Tatalina) 


CURVE  NO.  PLACE 

22  TATALINA 

23  TATALINA 

24  TATALINA 


REaECTIVITY  ( OB  ) 

FIGURE  32:  Comparative  Distributions  (Tatalina) 


REFLECTfVfTY  ( DB  } 

FIGURE  33:  Sea  Clutter  (Cape  Newenham) 


FIGURE  34:  Comparison  of  Rayleigh,  Log-Normal 

and  Contaminated-Normal  Distributions 


V.  SPECTRAL  DENSITY  MEASUREMENTS: 


Terrain  Clutter: 


Measured  spectral  densities  of  land  clutter  returns  are 
displayed  in  Photos  1 through  19.  Photos  1-3  give  a typical 
spectral  density  for  mountains  and  the  worst  measured  cases  for 
partially-wooded  hills  and  heavily-wooded  valleys.  For  use  in 
radar  design > an  attempt  was  made  to  describe  or  bound  the  worst 
case  two-sided  clutter  spectra  with  idealized  models  defined  by  a 
DC  impulse  and  a smooth  AC  function. 

P2Cf)  - Afig  (f)  + B GCf) 

For  mathematical  convenience,  the  DC  impulse  is  defined  as 

ScCf)  - I'  (I) 


where  ir  ^ { 


1>  l<  2 


0, 


fl>  L 
2 


c > normalizing  bandwidth 

In  determining  the  DC  impulse  component,  several  assumptions 
were  used.  First,  although  the  spectrum  analyzer  was  DC  coupled, 
it  did  not  display  DC.  Therefore,  all  land  clutter  measurements 
were  performed  with  the  peak  energy  of  each  return  adjusted  to 
within  3 dB  of  the  maximum  input.  Then,  the  energy  within  the 
filter  centered  at  zero  frequency  was  approximated  by 
2 


'DC 


( 


1 - 


If 


PRF 

2 lV(f)ldf  ) 


53 


Photo  1 


i'  1 0 0 H J 


'*  I L 'j 

-c  B 1 D B 


Photo  2 


Photo  3 


SPECTRAL  RETURN  FROM  MOUNTAINS  ( TYP  ) 


rt.  10  >>  4U  HLb  'w'Lb-lO 


SPECTRAL  RETURN  FROM  HILLS  ( WORST  CASE  ) 


SPECTRAL  RETURN  FROM  VALLEYS  ( WORST  CASE  ) 


Photo  5 


Photo  6 


Photo  7 


I 


5' 


HILLS 

Photo  14 


r 


I 

i 

j 


fe 


ti 


I 


i 


1 
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■where  A£ 


center  frequency  of  the  lowest  frequency  filter. 


Secondly,  the  spectrum  analyzer  internally  compensated  for 
spectra  amplitude  when  integration  times  were  changed.  Therefore, 
to  allow  a point-by-point  comparison  between  the  spectra  oscillo- 
graphs and  plot  ideal  functions,  the  normalizing  bandwidth  was  set 
at  e = .025  (the  narrowest  filter  of  the  spectrum  analyzer)  and 
the  total  energy  was  normalized  to 

P2(f)df  * .025 


Finally,  the  measured  spectral  spread  of  fixed  point  targets 
was  used  to  define  the  resolution  of  the  measurement  system. 

Since  this  resolution,  which  varied  from  .05  to  .25  Hz,  was  often 
much  larger  than  the  .025  Hz  resolution  possible  with  the  spectrum 
analyzer,  the  DC  coefficient,  A,  would  depend  on  the  choice  of  two 
sided  DC  filter  bandwidth  Bf.  Therefore,  Photos  1-3  were  graph! 
cally  analyzed  to  obtain  a DC-to-AC  ratio  a*  where 


V*DC 


Af  V*(f)df 


PRF 


V2Cf)df 


fl 


Bf 

2 


Figure  35  presents  a®  versus  filter  width  Bf  for  each  photograph. 
The  DC  coefficient  A can  now  be  given  as 


The  high  frequency  tails  were  compared  graphically  with  several 
standard  functions.  The  Gaussian  function,  suggested  by  several 
radar  publications,  did  not  adequately  describe  the  high  frequency 
tails.  Instead,  a more  descriptive  function  is  of  the  form^ 

""  ■ rr^ 

where  N - a constant 

ff.  * “3  dB  frequency. 

The  AC  coefficient  B can  be  given  as 

« .025 

B ■ 

(1  + a*)  E 

where  a*  ■ DC- to -AC  power  ratio 

- «0 

E - ' GCf)df 

•QO 

Descriptive  or  bounding  spectral  density  functions  for  moun- 
tains, partially-wooded  hills,  and  densely-wooded  valleys  are  given 
in  Figures  36  and  37.  The  typical  return  from  mountains  (peaks 
over  4,000  ft)  in  10  - 20  knot  winds  is  given  by  plot  1 (A).  The 
shape  of  the  measured  spectral  densities  from  mountains  was  uncer- 
tain due  primarily  to  the  unknown  instability  of  the  radar.  How- 
ever, the  measured  one-sided  spectral  density  PjCl  Hz)  was  always 
less  than  ‘ 50  dB  relative  to  the  maximum  energy  and  typically 
less  than  -60  dB.  A conservative  two-sided  spectral  density  which 
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boiinds  the  observed  values  for  mountains  is 


P2M  (f) 


0.964n 


where  Bf  ■ .067  Hz 


a»  - 27 


2.22  fc 

.2  Hz 


The  worst  case  spectra  measured  from  partially-wooded  hills 
(peaks  between  1,000  and  4,000  feet)  in  10  - 20  knot  winds  is 
given  by  plot  2 (X) . The  spectra  from  this  terrain  type  had  a 
narrow  peak  around  the  carrier  and  tails  that  decay  as  f“*  to  f* 
A conservative  two-sided  spectral  density  for  pax tially-wooded 
hills  is 


P2H  Cf) 


/ 0.00225  \ 


where  Bf  ■ . 2 Hz 

a*  ■ 9 

GCf)df  - 2.22 

«-00 

fc  - .5  Hz 

The  worst  case  spectra  measured  from  densely-wooded  valleys 
and  hillsides  (below  1,000  feet)  in  10  - 20  knot  winds  is  given  by 
plot  3 (0) , The  spectra  had  a broad  peak  about  the  carrier  with 
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significant  tails  decaying  as  f*  to  £'*  and  extending  beyond 
10  Hz.  A conservative  two-sided  spectral  density  of  the  form 


P2V  Cf) 


0.46  n 


.025 


0.00186 


*( 


(total  normalized  power  = .025) 
could  be  used  for  wooded  valleys  where 


Bf 

/ a* 

'CO 

GC£)df 


1 Hz 
. 85 

2.42  fc 

• 00 

£c  = 3 Hz 

B.  Sea  Clutter; 

Of  the  few  sea  clutter  spectra  recorded,  Photos  20  - 22 
present  the  worst  case  sea  clutter  spectra.  The  mean  doppler  is 
at  approximately  28  Hz  (6.4  knots)  and  the  3 dB  width  is  19  Kz  (4.4 
knots) . The  locally  measured  wind  conditions  were  10  knots  average 
at  the  11th  Weather  Squadron  Station  and  15-20  knots,  gusting  to 
25  knots,  at  the  radar  site.  Assuming  that  locally  measured  wind 
conditions  do  represent  the  wind  conditions  at  the  measured  range- 
azimuth  cell^,  a comparison  of  the  measured  spectra  with  values 
given  in  the  radar  literature  leads  to  the  following  observations: 

1.  The  spectral  shape  compares  favorably  with  the 

®It  is  recognized  that  local  wind  conditions  do  not,  in  general, 
describe  the  wind  conditions  several  miles  away. 


Photo  20 


Photo  21 


■«  ..  *■.  H ' I 


- f'l'M  ^ 

[ft  y “ I- 

V.V' 

• *»  v, 

I '» 


SEA  CLUTTER 


6<.l 


MEASURED  SPECTRA 


FIGURE  38:  Sea  Clutter  Power  Spectral  Density 


Gaussian  shape  often  suggested  (Figure  38) . 


2.  The  spectral  width  compares  favorably  with  other  experiment 
results  taken  under  similar  wind  conditions. 

3.  The  mean  doppler  was  1-2  knots  higher  than  other  pub- 
lished measurements  and  corresponded  to  wind  speeds  of  over 
40  knots.® 

4.  The  measurement  intervals  of  Photos  20  - 22  are  4,  10,  and 
34  seconds  respectively,  showing  little  change  in  spectral 
shape  over  these  integration  intervals. 

Beyond  approximately  20  nautical  miles  from  the  shore,  an  addi- 
tional spectral  component  was  noticed  (Photo  23).  Less  than  5 Hz 
wide  but  having  a maximum  doppler  of  at  least  70  Hz  (16  knots) , 
this  component  may  have  been  generated  by  large  ocean  swells  or 
gravity  waves. 

C.  Weather  Clutter; 


Measured  precipitation  doppler  spectra  are  given  in  Photos 
24  thru  37.  Most  of  the  precipitation  doppler  spectra  recorded  in 
the  clutter  study  were  taken  at  Cold  Bay  during  a gale.  Surface 
observations  taken  during  the  measurements  indicated  rainfall  rates 
of  up  to  8 mm/hour  at  Cold  Bay  and  over  28  mm/hour  at  Fort  Yukon. 
Radiosonde  data  taken  within  48  hours  of  the  measurement  period  at 
Cold  Bay  indicated  cloud  decks  (relative  humidity  >601)  of  10,000 
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to  16,000  ft  and  average  wind  speeds  of  up  to  42  knots  within  the 
cloud  deck  (see  Table  IV] . Surface  measurements  indicated  average 
wind  speeds  up  to  36  knots  with  gusts  to  44  knots. 


TABLE  IV: 

Measured  Wind  Velocities  and 
For  Cold  Bay  Measurements 

Directions 

9 HOURS  BEFORE  MEASUREMENTS 

HEIGHT  (FT) 

DEGREES 

KNOTS 

3700 

165 

42 

8777 

190 

28 

17154 

235 

36 

22435 

230 

46 

3 HOURS  AFTER  MEASUREMENTS 

HEIGHT  (FT) 

DEGREES 

KNOTS 

3402 

ZOO 

18 

8396 

205 

30 

16601 

205 

18 

217S3 

215 

30 

Photos  24  thru  30  were  taken  at  Fort  Yukon  during  rainy  and 
windy  conditions.  Photos  28  - 30  were  taken  within  5 minutes 
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of  each  other  and  show  the  relative  backscatter  of  rain  versus 
polarization.  Rain  backscatter  using  vertical  polarization  and 
circular  polarization  on  transmit  and  receive  were  less  than  back- 
scatter using  horizontal  polarization  by  approximately  8 and  13 
dB  respectively.  (Observation  of  horizontal  polarization  of  sev- 
eral minutes  before  and  after  this  series  was  taken  revealed  little 
amplitude  or  frequency  change  of  the  spectra.  Therefore,  the  differ- 
ence in  the  three  spectra  is  assumed  to  be  due  to  polarization.) 

Photos  31  - 37  were  taken  at  Cold  Baya.  Photos  36  and  37  give 
a good  example  of  the  spectral  environment  within  30  nautical  miles 
of  Cold  Bay  during  the  measurement  period.  Photo  36  is  an  averaged 
spectrum,  taken  over  3.6  second  interval^,  of  the  backscatter  from 
a range-azimuth  cell  over  land  at  a range  of  25  NM.  Photo  37  is  an 
instantaneous  spectrum,  taken  over  a 2 second  interval^,  of  the 
backscatter  from  a range-azimuth  cell  over  the  Bering  Sea  at  a range 
of  20  NM. 

^In  reading  the  Cold  Bay  precipitation  spectra  photos,  note  that 
the  radar  pulse  repetition  frequency  (PRF)  was  244  Hz  causing 
spectral  foldover  at  122  Hz. 

^For  "instantaneous"  spectra,  the  analyzer  requires  a time  interval 
Ti  which  can  be  given  by  Ti(secs)  = 1024  samples/spectrum  ^ (2,56)fD 
where  fj)  = maximum  display  frequency. 

For  "averaged"  spectra,  the  time  interval  can  be  given  by 
T^Csecs)  = Ti  X [ 1 + ] 

where  N = number  of  spectra  to  be  averaged. 

Note  that  successive  501  redundant  spectra  are  averaged. 


Table  V gives  a comparison  of  several  photos  with  values 
given  by  Nathanson^.  Since  some  pictures  were  not  symmetrical 
or  Gaussian,  a BW.3JB  and  a BW-iq^b  is  given.  Despite  the  long 
averaging  time,  a Gaussian  function  centered  at  the  measured 
median  values  with  a maximum  equal  to  the  mode  values  and  having 
the  standard  deviations  given  by  Nathanson  bound  all  measured 
values . 
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D.  False  Alarms: 

A table  of  false  alarms  for  each  site  is  given  in 
Table  VI.  The  three-pulse  MTI  canceller  of  the  FPS-93A  was 
calibrated  with  respect  to  velocity  shape  factor  and  cancellation 
by  the  ADC/AAC  team.  Then,  using  a PPI,  the  false  alarms  thru  the 
MTI  channel  (without  STC)  were  studied.  Nearly  all  the  studied 
false  alarms  were  tracked  several  scans  and  moved  measurable 
distances,  thereby  having  the  characteristics  of  targets  rather 
than  noise.  These  false  alarms  have  been  defined  as  angels.  The 
number  of  false  alarms  were  obtained  by  counting  the  number  of 
false  alarms  in  an  average  populated  sector  of  45“  and  multiplying 
by  8. 


^Nathanson,  F.  E.,  Radar  Design  Principles.  McGraw-Hill  Book  Co., 
Inc.,  New  York,  NY,  1969. 
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TABLE  V:  Measured  Precipitation  Spectra 


Photo 
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Photos  38  thru  55  show  the  spectra  of  the  two  angel  categories 
measured.  The  angels  in  one  category  are  called  dot  or  point 
angels  because  they  occupy  only  one  range  cell  and  resemble  small 
targets. 

The  dot  angels  were  observed  to  have  the  following  character- 
istics : 

a.  Direct  measurement  of  the  amplitudes  of  the  dot  angels 
was  not  possible  due  to  ground  clutter.  The  attenua- 
tion required  to  clear  the  CTI  channel  PPI  display  of 
all  dot  angels  is  given  in  Table  VI.  Analysis  of 
spectral  photos  of  some  of  the  returns  reveal  signal- 
to-noise  ratios  up  to  30  dB  in  a 200  KHz  video  band- 
width. 

b.  Distributions  of  measured  doppler  velocity  typically 
ranged  between  ^ 25  knots.  The  highest  doppler 
measured  was  42  knots.  The  spectrum  of  dot  angels 
always  resolved  to  one  or  more  approximations  to  line 
spectrum  at  some  time  during  the  period  of  observation. 

Attempts  to  associate  the  dot  angels  with  visual  observations 
of  targets  were  unsuccessful.  However, areas  with  many  radar 
angels  per  scan  were  also  areas  with  large  bird  populations.  There- 
fore, it  would  not  be  unreasonable  to  assume  the  angel  cross-sections 
to  be  that  of  birds,  such  as  ducks,  geese,  hawks,  and  sea  gulls. 
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The  other  category  of  angels  is  called  distributed  angels 
because  they  usually  occupied  more  than  one  range  cell.  These 
angels  were  identified  by  the  following  characteristics: 

(1)  They  occur  night  or  day. 

(2)  They  do  not  seem  to  move  radially  or  azimuthally  over 
periods  of  several  minutes  but  have  wide  noise- like 
doppler  spectra. 

(3)  The  signal  strength  fluctuated  several  dB  over  a period 
of  minutes  but  usually  was  very  weak.  Two  exceptions 
are  shown  in  Photos  48  thru  55. 

(4)  They  occurred  during  clear  weather. 

(5)  No  obvious  exj)lanation-waa  evident.  However,  clear  air 
turbulence,  aurora,  and  covies  of  birds  or  insects  are 
possible  candidates. 

(6)  Less  than  1%  of  the  angels  observed  were  classified  under 
the  category  of  distributed  angels. 

Other  false  alarms  noted  during  the  study  were  those  due  to 
anomalous  propagation.  These  returns  were  correlated  with  moun- 
tains or  islands  up  to  150  NM  away  but  had  post-MTI  signal-to- 
noise  ratios  of  over  30  dB,  Unfortunately,  the  anomalous  propa- 
gation conditions  did  not  occur  during  the  spectral  analysis 
intervals  and  therefore  no  spectral  data  could  be  taken. 
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Photo  41 
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DISTRIBUTED  ANGELS 
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VI.  SUMMARY  AND  CONCLUSIONS: 


Clutter  measurements  of  the  Alashan  radar  environment  have 
been  presented  in  terms  of  several  amplitude  cumulative  distribu- 
tions of  land  and  sea  clutter  effective  reflectivity  coefficient 
Co'q3  and  the  power  spectral  density  of  land,  sea,  weather  and 
angel  clutter.  Models  have  been  given  for  L-Band  land  clutter 
for  use  in  radar  design. 

The  discussion  on  land  clutter  data  presents  the  following 
conclusions : 

1.  Homogeneous  clutter  whose  scatterers  are  larger  than  a 
carrier  wavelength  and  much  less  than  a resolution  cell  has  log- 
normally  distributed  reflectivity  coefficient.  The  smaller  the 
grazing  angle  the  larger  the  spread  between  the  50  and  the  84 
percentile, 

2.  The  more  heterogeneous  the  data  becomes,  where  large 
areas  within  a sector  become  shadowed  and  large  exposed  areas 
become  prominent,  the  distribution  becomes  more  like  a contami- 
nated-normal with  very  irregular  knees  and  bends. 

3.  As  the  areas  of  shadowing  and  exposure  become  numerous 
by  observing  larger  and  larger  sectors,  the  distribution  takes  on 
a smooth  transition  from  a small  slope  at  the  50  percentile  to  a 
large  slope  at  the  low  percentiles. 
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4.  The  tails  of  all  distributions  approach  a log-normal 
function  at  the  low  percentiles. 

Measurements  of  land  clutter  spectra  have  shown  that  the 
Gaussian  function  often  used  in  radar  design  does  not  'idequately 
describe  the  high  frequency  tails  of  the  measured  spectral  dis- 
tribution function.  A more  accurate  description  is  provided  by 
the  functions 
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where  N is  typically  between  3 and  S and  f(.  represents  the  -3  dB 
point  of  the  AC  curve.  The  AC  -3  dB  spectral  widths  f^  and  DC-to- 
AC  ratios  (a*)  of  terrain  clutter  returns  ranged  from  fc  < .25  Hz 
and  a*  » 27  for  mountains  to  fc  “ 3 Hz  and  a*  < 1 for  valleys. 

The  worst  case  measured  sea  clutter  spectra  displayed  nearly 
Gaussian  shaped  spectral  widths  close  to  what  is  presented  in  the 
literature,  but  much  higher  mean  doppler  than  locally  measured  - 
wind  conditions  would  have  inferred.  Also,  a high  mean  doppler, 
narrow  spectral  component  was  observed  beyond  20  NM  and  is  presumed 
to  be  from  ocean  swells. 

The  measured  weather  spectra  displayed  means  of  up  to  at 
least  the  unambiguous  velocity  of  the  measurement  radar  (28  knots) 
and  various  spectral  shapes  bounded  by  the  design  values  given  in 
standard  radar  texts. 
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Finally,  the  false  alarms  occurring  through  the  MTI  channel 
at  each  site  were  presented.  Most  of  the  false  alarms  not  caused 
by  weather,  sea  clutter,  land  clutter,  or  noise,  were  termed  angels 
and  were  attributed  to  birds. 
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APPENDIX  A 

DERIVATION  OF  REFERENCE  THRESHOLDS  FOR  ac  AND  o^: 

The  reference  thresholds  for  the  effective  radar  cross- 
section  0^  and  the  effective  normalized  radar  cross-section  o© 
were  derived  from  the  standard  radar  range  equation  and  the  para- 
meters of  the  calibrated  threshold  comparator  used. 


The  effective  radar  cross-section  can  be  defined  as: 
PrR'  I ^ I (4Tr)* 

[PtOtGrl*  J 

[^]  [’■]  [ssfc][(r)' ■■■’■] 

[PtG*G_X*  “I 
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+ 10  Log  ^1)  (4Tr)’  j + LattN  " ^rc 
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10  Log 

[vi) 
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®rdBm 

®cdBm*  “ ^ 


rdBm 


30  dB 


+ 40  Log  T - 10  Log 
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CA-l) 


(A- 2) 


+ 330  - Grc  + Latin 

The  threshold  varies  as  a function  of  time  such  that 
Sr  dBm  + 40  Log  T « Constant  = K 


return  from  a constant  radar  cross-section  as  a 
function  of  range. 


Let  reference  T^q  and  S;f(10)  be  given 
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Therefore,  the  reference  threshold  is  given  as 
168.7  - 10  Log  [ 
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Similarly,  the  reference  threshold  for  the  effective 
normalized  radar  cross-section  can  be  defined  as 
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The  threshold  varies  as  a function  of  time  such  that 
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return  from  a constant  reflectivity  as  a function  of 
range. 


Let  references  and  Sr(lO)  be  given  as 
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Therefore,  th^  reference  threshold  is  given  as 
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APPENDIX  B 


I.  ERROR  ANALYSIS  OF  MEASUREMENTS; 

A.  Amplitude  Measurements: 

The  accuracy  of  the  method  described  depends  upon  the 
accuracy  to  which  the  radar  parameters  can  be  measured  or  esti- 
mated and  the  extent  to  which  the  assumption  of  the  receiver 
linearity  holds  true®. 


The  assumed  accuracy  of  the  radar  parameter  measurements 
are  given  in  Table  B-I.  The  receiver  gain  error  is  given  as  .5 
dB  rms.  This  was  found  to  be  a good  approximation  up  to  "he  1 
dB  compression  point  of  the  receivers  at  each  site.  The  dynamic 
range  of  the  receivers  at  three  of  the  sites  are  given  in  Figure 
B-1.  x’he  other  two  sites  had  modifications  which  allowed  attenu- 
ation to  be  placed  prior  to  the  first  active  device,  eliminating 
any  dynamic  range  problems.  An  approximate  relationship  between 
the  1 dB  compression  point,  Oq  max,  and  range  was  used  to  deter- 
mine the  areas  where  valid  data  could  be  obtained^. 

®Due  to  time  and  operational  restrictions,  the  use  of  a calibrated 
target  was  not  feasible.  Therefore,  the  radar  range  equation  was 
assumed  to  exactly  describe  the  measurement  problem. 

^Despite  the  care  taken  at  each  site,  a few  measurements  had 
returns  exceeding  the  one  (1)  dB  compression  point  of  the 
receiver.  Listed  in  Tables  B-IV  - B-VII  are  the  calculated 
maximum  values  of  effective  radar  cross-section  or  normalized 
radar  cross-section  (omax)  which  the  error  analysis  applies. 
Also  listed  is  the  region  in  the  tails  of  the  plotted  distribu- 
tion C%  min)  corresponding  to  this  Omax-  Beyond  this  value, 
the  error  increases  due  to  compression  and  the  tails  become 
distorted. 


TABLE  B-I 


Accuracy  of  Radar  Parameter  Measurements 


PARAMETER 

RMS 

ERROR 

RMS 

ERROR 

C.dBl  .CBJ. 

Average  Power  (Pavgl 

— 

.5  dB 

Pulse  Width  (t) 

.05 

.22  dB 

Pulse  Repetition  Rate 

.5  Hz 

.01  dB  max 

Wavelength  (A) 

.OOlX 

.004  dB 

Transmitter  RF  Loss  (L^) 

— 

.5  dB 

Receiver  RF  Loss  (Lr) 

— 

.5  dB 

Propagation  Loss  (Lp) 

— 

Variable 

Radome  Loss  CLra) 

— 

.5  dB 

Sine  3 

.001 

.18  dB 

Threshold  Constants 

a.  R^  Threshold  (K) 

— — n 

Variable 

b.  R^  Threshold  (C) 

— 

Variable 

Receiver  Gain  (Grc) 

— 

.5  dB 

Attenuator  Pad  Loss  (LatTN) 

— 

.1  dB 

Antenna  Gain  (Gt»  Gr) 

— 

Variable 
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ERROR 

DISTRIBUTION 
Log  Normal 


Log  Normal 
Log  Normal 

Log  Normal 


Log  Normal 
Log  Normal 
Log  Normal 


FIGURE  B-1:  L-Band  Mixer  - Preamplifier  Dynamic  Range 


The  antenna  electrical  tilt  angle  at  each  site  could  not  be 
changed  due  to  operational  considerations.  Therefore,  a curve  of 
the  measured  antenna  gain  function  in  the  vertical  plane  provided 
by  the  antenna  manufacturer  had  to  be  used  CPiEure  B-2).  For  the 
purposes  of  this  exercise.  Figure  8-2  was  assumed  to  exactly 
describe  the  antenna  at  each  of  the  sites.  From  the  vertical 
cuts  of  the  horizontal  pattern,  the  azimuthal  beamwidth  was 
found  to  be  1,35*  ^ .05.  Therefore,  Sin  ©p  ■ .024  ^ .001.  If 
the  electrical  tilt  angle  of  the  antenna  and  the  depression  angle 
of  the  ray  directed  toward  a clutter  patch  were  known,  the  antenna 
gain  could  be  read  off  the  curve.  The  tilt  angles  were  measured 
within  the  tolerance  given  in  Tables  B-IV  thru  B-VIII.  The 
depression  angle  toward  a clutter  area  can  be  found  by 


Sin"  * 


- ( h2  - hi)_  ( hz  - hi)  ^ R 1 


where  h2 
111 
R 
re 


height  of  clutter  area  (ft) 
height  of  radar  (ft) 
range  to  clutter  area  (ft) 
true  radius  of  earth  (ft) 


The  ^ earth- radius  approximation  to  the  effects  of  standard  re- 
fraction was  assumed  in  the  derivation  of  the  above  equation. 
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FIGURE  B- 


: FPS-93  Antenna  Pattern  Principal 
Plane  Elevation  Horizontal  Polari- 
zation 1300  MHz 
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Since  most  areas  measured  have  a range  and  height  extent, 
minimum  and  maximum  depression  angles  were  determined  for  each 
area.  After  conversion  to  gain  through  the  use  of  Figure  B-2, 
an  average  gain  (G) , maximum  error  (Epk) , and  RMS  error  (Erms) 
(assuming  log  normal  distribution)  was  associated  with  each 
(ErMS)  measurement.  For  convenience,  Epk  £ 2o  ■ 2 Erms.  The 
total  antenna  RMS  gain  error  was  determined  by 

Eg  ■ [ E*rms  + g*  (Tilt  Angle) 

where  g*  (Tilt  Angle)  - Squared  Error  in  Gain  due  to  Tilt  Error 

« [3.8  (Tilt  Error)  ]* 

The  propagation  loss  error  is  determined  by  the  range  extent 
of  the  clutter  area  being  measured  and  can  be  found  from  Figure 
B-3. 


The  threshold  error  was  defined  as  the  RMS  deviation  from  the 
desired  constant.  The  measured  error  versus  range  is  given  in 
Tables  B-II  and  B-III  for  the  oq  threshold  and  Oq  threshold, 
respectively.  As  one  can  see,  the  mean  error,  which  was  included 
into  the  calculation  of  reference  parameter  for  each  measurement, 
and  the  standard  deviation  of  the  error  about  the  mean  depend  on 
the  range  and  range  extent  of  the  clutter  area.  These  are  listed 
for  each  measurement  in  Tables  B-IV  thru  B*VIII. 
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Since  it  is  unlikely  that  all  the  errors  vrill  add  linearly, 
an  RMS  value  was  considered  more  representative.  The  RMS  error 
of  the  absolute  value  for  the  oc  and  oo  calculations  can  be  de- 


rived from  equations  (A- 3)  and  (A- 7),  respectively. 


Ek  ♦ + Ex  * ♦ 4Eq  + 4Ex  + E^p  + EGyc 


* ^LattN  + ELx  + Elra 


Ec  + ®PavG  * ^ 4Ex  + Ej^p  + Bsineg 


* ^^ATTN  * ^^RC  * ^Lj 


where  E ■ standard  deviation  in  dB 


Pt  - PaVG  - ■fdBsgc 
Ls  - EpLpL^Lpy^ 


PWdBHs 


Substituting  the  values  given  in  Table  B-I,  the  mean  squared 


errors  are 


Ep-  - 1.3  + Ejc  + 4Eg  + Elj 


l _2 


Eg'  - 1.3  + Ej.  + 4Eg  + Et 


These  values  for  each  measurement  are  given  in  Tables  B-IV  thru 
B-VIII. 


TABLE  B-II 


ac  Threshold  Error 


Sr 

Teal 

Tteeas 

S TdBsec 

3 K 

(dBn^ 

Cusec) 

Cusec) 

(dB) 

(dB) 

- 5.0 

123.5 

145,7 

.72 

2.87 

-12.0 

185.25 

197.0 

.27 

1.07 

-17.0 

247.0 

247.0 

0 

0 

-24.1 

370.5 

361.5 

-.11 

-.43 

-29.1 

494.0 

476.1 

-.16 

-.64 

-33.0 

617.5 

596.6 

-.15 

-.60 

-36.1 

741.0 

716.1 

-.15 

-.59 

-38.8 

864.5 

835.6 

-.15 

-.59 

-41.1 

998.0 

958.1 

-.13 

-.53 

-43.2 

1111.5 

1094.0 

o 

• 

1 

00 

CM 

1 

-45.0 

1235.0 

1229.0 

CM 

o 

1 

1 

O 

00 

-52.0 

1852.5 

1906.0 

+.12 

+ .50 

-57.0 

2470.0 

2581.0 

+.19 

+.76 

I ‘ 


TABLE  B-III 


Oo  Threshold  Error 


Sr 

T cal 

T neas 

a TdBsec 

d C 

(dBm) 

(usee) 

(usee) 

(dB) 

(dB) 

-15.0 

123.5 

144.5 

+ .68 

2.05 

-20.3 

185.25 

20C.8 

+.35 

1.05 

-24.0 

247.0 

241.7 

-.09 

-.28 

-29.3 

370.5 

354.1 

-.0 

-.59 

-33.1 

494.0 

464.1 

-.27 

-.81 

-36.0 

617.5 

581.7 

-.26 

-.78 

-38.3 

741.0 

701.1 

-.24 

-.72 

-40.4 

864.5 

824.6 

-.21 

-.62 

-42.1 

988.0 

936.3 

-.23 

-.70 

-43.6 

1111.5 

1074.0 

-.15 

-.45 

-45.0 

1235.0 

1221.0 

-.05 

-.15 

-50.3 

1852.5 

1938.0 

+ .20 

+ .59 

-54.0 

2470.0 

2549.0 

+.14 

+.41 
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TABLE  B-IV:  Measurement  Errors  - Fort  Yukon 


♦Error  in  Tilt  Angle  Measurement 


Measurement  Errors  - Tatallna 


•Error  in  Tilt  Angle  Measurement 


Measurement  Errors  - Tatalina  (Continued) 


oo  ooooooo 


« m tM  o» 

+ + I 


^0  M K) 


o lO 

»-•  rH  rH 

♦ + + 


r in 

CM 

CM 

4 

pH 
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o (-imiAi-it-iiotAr-itniOi-imLni-tcM 

O 
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to 
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to 

to 

o 

• 

CM 
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CM 

CM 

• 

• 

• 

• 

• 
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•Error  in  Tilt  Angle  Measurement 


TABLE  B-VIII:  Measurement  Errors  - Cold  Bay 


TABLE  B-VIII:  Measurement  Errors  - Cold  Bay  (Continued) 


•Error  In  Tilt  Angle  Measurement 


B.  Spectral  Measurements : 


The  errors  in  the  spectral  measurements  were  essentially 
the  spectral  contamination  contributed  by  the  interface  equipment 
and  the  radar. 

The  spectral  contamination  from  the  interface  equipment 
(excluding  £he  30  MHz  local  oscillator)  was  essentially  the  ampli- 
tude distortion  from  the  ~ weighting  of  the  sample- and-hold 
circuitry  and  a low  pass  filter  as  shown  previously  in  Figure  9. 
Distortion  measurements  were  made  to  determine  the  regions  of 
linearity  in  the  video  circuitry  using  two  equal  amplitude  audio 
tones.  Harmonic  distortion  contributed  by  the  spectrum  analyzer 
was  immeasurable  (better  than  -74  dB  below  the  amplitude  of  one 
tone)  for  signal  levels  below  A/D  converter  saturation.  The  max- 
imum 2nd  harmonic  distortion  within  0 - 300  Hz  contributed  by  the 
200  KHz  amplifier  and  sample-and-hold  was  better  than  -48  dB  below 
the  fundamental  for  signal  inputs  of  250  mV  pk-pk  or  less.  Third 
order  distortion  for  this  input  was  better  than  -55  dB  for  fre- 
quencies below  100  KHz.  Most  measurements  were  performed  with 
inputs  less  than  250  mV  pk  - pk.  The  linear  limit  amplifier  and 
heterodyne  detector  contributed  inter-modulation  distortion  less 
than  -50  dB  for  video  outputs  of  250  mV  pk  - pk  or  less.  At  sites 
where  attenuation  could  not  be  placed  ahead  of  the  mixer-preamp, 
only  signals  at  least  IS  dB  below  the  1 dB  compression  point  were 
analyzed  . 
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It  was  assumed  that  signals  at  this  level  or  below  would  contri- 
bute negligible  spectral  distortion. 

The  spectral  contamination  caused  by  the  transmitter  consisted 
primarily  of  spurious  frequencies  and  spreading  of  observed  spectra 
due  to  oscillator  instability.  A spectral  analysis  of  the  FPS-65 
transmitter  at  Verona  Test  Site  (closed  loop)  disclosed  a random 
noise  component  and  discrete  frequencies.  The  random  noise  com- 
ponent was  centered  about  the  carrier  and  was  at  least  60  dB  down 
from  the  carrier  at  .1  Hz  and  could  be  ignored.  The  discrete 
frequencies  with  amplitudes  above  -60  dB  were  60  Hz  (-50  dB)  and 
a spurious  (-58  dB)  which  was  caused  by  intermodulation  of  the 
PRF  with  60  Hz.  When  observing  nearby  "stable"  clutter  returns, 
the  spurious  frequencies  modulated  each  return  providing  a method 
for  identifying  this  form  of  contamination.  For  convenience,  this 
was  the  method  used  during  the  Alaskan  measurements  to  determine 
transmitter  spurious  frequencies. 

The  short-term  stability  of  the  radar  oscillators  and  timing 
system  was  the  limiting  factor  determining  the  resolution  of  the 
spectral  analysis.  The  changes  in  oscillator  phase  between  trans- 
mission and  reception  from  clutter  often  resulted  in  a spectral 
spreading  or  smearing  of  several  tenths  of  a Hertz.  For  instance, 
a .25  Hz  shift  was  imposed  on  a return  300  usee  away  at  one  site 
due  presumably  to  the  stalo  oven  cycling.  In  many  cases,  however, 
the  spectral  spread  was  less  than  .05  Hz.  Therefore,  the  resolu- 
tion of  the  measurements  varied  but  was  less  than  1 Hz  in  ail  cases. 
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The  spectrum  analyzer  computes  and  displays  the  one-sided 
energy  spectral  density  of  the  input  waveform.  However,  it  does 
not  display  the  zero  frequency  filter.  Therefore,  an  indirect 
method  was  used  to  determine  the  DC-to-AC  power  ratio  of  land 
clutter.  The  spectrum  analyzer  is  peak  voltage  limited  and 
internally  relates  all  spectral  amplitudes  to  this  maximum  level. 

By  setting  the  input  signal  to  just  below  saturation  when  the  coho 
was  in-phase  with  the  return,  the  scale  was  calibrated  to  the  max- 
imum peak  energy  of  the  heterodyned  return.  Then,  the  relative 
energy  in  the  zero  filter  was  approximated  by 

Vdc  = C 1 ’ ^ lVCf)!df  )* 

This  approximation  worked  well  when  the  DC  energy  is  much  greater 
than  the  AC  energy  and  the  selected  filter  bandwidth  Bf  was  larger 
than  the  spectral  spread  due  to  instability. 

Finally,  the  term  actually  measured  by  the  equipment  was 
average  voltage  or  average  energy  spectral  density  over  a specific 
time  interval.  The  measured  process  would  have  to  be  stationary  to 
have  a correlation  function  and,  consequently,  a power  spectrum. 
Although  radar  clutter  is,  in  general,  a non- stationary  process, 
it  can  usually  be  assumed  to  be  stationary  over  the  observation 
interval  of  a long  range  search  radar.  This  allows  a stringent 
stationary  clutter  model  to  be  constructed  for  the  purpose  of 
comparing  different  radar  designs.  Since  the  purpose  of  this 
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study  was  to  obtain  data  for  such  a clutter  model,  the  term  power 
spectral  density  is  used  throughout  the  report  to  describe  the 
measured  data  as  well  as  the  idealized  functions. 


APPENDIX  C 


AVEt^ ^GE  Op  AND  DEPRESSION  ANGLES  FOR  LAND  CLUTTER  MEASUREMENTS: 

Table  C>I  lists  the  average  depression  angle  and  the  average 
effective  reflectivity  coefficient  Oq  for  the  measurement  areas 
presented  in  Section  IV.  Unless  otherwise  noted,  the  6 usee 
measurements  were  used  in  the  calculations. 

In  order  to  obtain  a histogram  approximation  to  the 
density  function  was  obtained  from  the  original  data  of  each 
measurement.  All  the  range  cells  in  the  measurement  area,  in- 
cluding those  range  cells  with  immeasurable  returns  due  to  the 
presence  of  noise,  were  included  in  the  sample  size  to  obtain 
the  true  average  value.  Since  the  actual  distribution  of  the 
immeasurable  returns  is  unknown,  two  extreme  distributions  were 
considered.  In  the  first  case,  the  immeasurable  returns  were 
assigned  a Oq  just  below  the  threshold  sensitivity  of  the  measure- 
ment equipment  with  a probability  of  occurrence  equal  to  the 
percentage  of  range  cells  with  an  immeasurable  return.  In  the 
second  case,  the  immeasurable  returns  were  assigned  a “ 0. 

These  cases  provide  an  upper  and  lower  bound  to  the  true  average 
value. 

With  one  exception,  the  average  values  calculated  in  both 
extreme  cases  differed  by  less  than  .1  dB.  The  exception  was 
the  measurement  of  the  Fort  Yukon  swamp  where  a difference  of 
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.5  dB  was  calculated.  Compared  with  the  accuracies  of  the 
measurements,  these  differences  are  negligible.  Therefore, 
the  actual  distribution  of  the  immeasurable  returns  is 
irrelevant. 


TABLE  C-I:  Average  Reflectivity  and  Depression  Angles 


SITE 

TERRAIN 

TYPE 

AREA 

CODE 

AVERAGE 
DEPRESSION 
ANGLE  (Deg) 

(dB) 

(t  “ 6 usee) 

Fort  Yukon 

Swamp,  Muskeg 

1 

.15 

-60 

High  Hills,  Mtns. 

3 

0 

-37 

Hills 

4 

0 

-24 

Hills 

6 

0 

-17 

Indian  Mt. 

Hills 

2 

.6 

-22 

Hills 

3 

.6 

-14 

Hills 

4 

.4 

-25 

Hills 

5 

.7 

-16 

Mtns . 

6 

.6 

-28 

Mtns . 

7 

.6 

-210 

Tatalina 

Everything 

2 

.7 

-20 

Everything 

3 

.7 

-23 

Hills 

6 

.6 

-19 

Muskeg , Swamp 

7 

1.2 

-17 

Mtns . 

8 

.1 

-18 

High  Hills,  Mtns. 

9 

.3 

-19 

Cape  Newenham 

Hills 

3 

.5 

-41 

Hills 

5 

.6 

-32 

Cold  Bay 

Volcanic  Mtns. 

1 

- .8* 

-26 

Volcanic  Mtns. 

2 

- .8* 

-30 

(+)  -3  usee  pulse  width 

* -negative  depression  angle  indicates  angle  above  tangent 
to  surface 


